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Abstract 

We study “field space entanglement” in certain quantum field theories consisting of N number 
of free scalar fields interacting with each other via kinetic mixing terms. We present exact analytic 
expressions for entanglement and Renyi entropies between arbitrary numbers of scalar fields by 
which we could explore certain entanglement inequalities. Other entanglement measures such 
as mutual information and entanglement negativity have also been studied. We also give some 
comments about possible holographic realizations of such models. 
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1 Introduction 

Quantum entanglement offers different measures to capture some non-local properties in quantum 
field theories (QFTs). There are various measures for quantum entanglement including entangle¬ 
ment and Renyi entropies [1] which measure the amount of quantum entanglement between various 
parts of the Hilbert space of the theory. Among these measures, specifically entanglement entropy 
(EE) has recently gained a huge amount of interest. 

In this context, the most common way available in the literature for studying quantum en¬ 
tanglement is based on a one-to-one correspondence between localized degrees of freedom of local 
quantum field theories and plane waves as a particular complete basis spanning their total Hilbert 
space. Based on such a map the Hilbert space is decomposed as 71 = 71 a ^7-Lb, where A and B 
correspond to spatial subregions such that A4 = AuR is a constant time slice of the manifold which 
the QFT is defined on. Such a decomposition is reliable up to the spatial resolution introduced by 
the UV cut-off of the theory. The spatial subregions A and B are defined via a co-dimension-two 
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surface dA. Following such a decomposition and tracing out either part A ox B leads to a measure 
for the quantum entanglement between localized degrees of freedom in spatial regions A and B. 
We denote this type of EE as “spatial entanglement entropy” (SEE). Some well known features of 
entanglement entropy such as the celebrated area-law divergence [2,3] is peculiar to SEE. 

SEE is not the only type of EE one can define between various degrees of freedom of a single field. 
There are other types of EE corresponding to different Hilbert space decompositions. For example 
one can decompose a given Hilbert space into states with specific energies and consider the EE 
referring to given scale of energy A. This type of EE is known as “momentum space entanglement 
entropy” which measures the EE between degrees of freedom of a single field below and above a 
given energy scale A in the momentum space (see e.g. [4]).^ 

If more than one field lives in a field theory, one may ask about probable entanglement between 
degrees of freedom corresponding to different fields. In contrast to various EE measures defined 
between different degrees of freedom of a single field, the entanglement between degrees of freedom 
of different helds is caused via possible interactions between them.^ Using the terminology of 
reference [6], we denote this type of EE as “held space entanglement entropy” (FSEE). 

It is worth to note that Ryu-Takayanagi proposal [9-12] for holographic entanglement entropy 
is by construction a proposal to compute SEE in a held theory which supports classical Einstein 
theory as a gravity dual. A natural question which may arise is about the possibility of a holographic 
realization for other types of EE e.g. FSEE. We are not going to answer this question in this paper 
and we will only give some comments about it in the section 5. Recently some arguments about 
this interesting question has appeared in the literature specihcally in [13] and [6] (see also [14] for 
some related holographic improvements). 

In this paper we try to further investigate the notion of ESEE from a held theoretic point of view. 
To do so we consider various held theories which are interacting with each other. The interaction 
between these held theories is responsible for generating entanglement between them. In order to 
study the entanglement between these theories we integrate out a generic number of them which 
leads to a reduced density matrix. Next we follow the standard procedure to study entanglement 
and Renyi entropies. 

Eor simplicity we focus on scalar held theories with Gaussian interactions between them. Since 
such models are Gaussian, they are analytically tractable to a satisfactory extent, and thus we 
consider them as a simple laboratory to study some general properties of ESEE. Explicitly we work 
out the generic reduced density matrix of such models and study entanglement and also all Renyi 
entropies analytically. A similar construction have been previously studied in [13] and in the context 

^ There are also two other types of entanglement discussed in the literature: The first one which is called “entan¬ 
glement in theory space” is defined via gauging (un-gauging) two theories with global symmetries in [5]. We would 
like to thank Mukund Rangamani for bringing our attention to this reference. The other one which is called “global 
symmetry entanglement” is defined via partitioning the symmetry group in [6]. 

^We are aware of some studies which can be considered as quantum mechanical counterparts of such an analysis, 
including reference [7] where entanglement between non-interacting qubits is studied and also reference [8] where a 
particle partitioning is considered for studying entanglement entropy. 
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of condensed matter physics in [15-18]. 

The authors of reference [13] have considered two free scalar field theories denoted by cj) and iIj 
which interact homogeneously in a d-dimensional space-time via two types of interactions; kinetic 
mixing (marginal) and massive interactions. They have decomposed the total Hilbert space of the 
theory as Ti = (g) and integrated out the states in and worked out the entanglement and 

Renyi entropies between (j) and ^ in the ground state which is no more a direct product due to the 
interaction between them. 

In this paper we generalize the procedure of reference [13] in the sense that we consider N free 
field theories defined on a common d-dimensional flat space-time which interact with each other. 
The action is thus given by 

^ ~ J '^2 (02) + ■ ■ ■ + Ttv (0Ar) + Tint, (0i)] ; (1-1) 

where T* (0i) with i = 1, 2, • • • , iV denote the Lagrangian density of free held theories and Tint. {4>i) 
denotes all possible interactions between them. We are interested in entanglement and Renyi 
entropies between these held theories which is generated via the interaction term Tint, (0i)- The 
total Hilbert space of this model can be decomposed as 


"^tot. = ® ^2 <8) • • • (8) Ti-N, 


where T-LiS denote the Hilbert space of each held theory dehned by Tj (0*). We are interested in the 
entanglement between generic m number of these held theories with the rest {N — m) of them. To 
do so we consider the following more compact notation for the decomposition of the total Hilbert 
space as 

T-Ltot. = T-i(^rn) ® T~{-{N-m) (1-2) 

where ^(m) is dehned as = 'Hi 8> H 2 ® • • • 8) Hm and H(Ar_m) similarly denotes the Hilbert 

space of the rest {N — m) held theories. In such a way we dehne the reduced density matrix 
by tracing out the H(jv-m) part of the Hilbert space 


P(m) [Ptot.] ) 


which leads to the following dehnition of entanglement and Renyi entropies 


Sent, (m) = -Tr [p(^) log p^rn)] , (m) = log Tr 


P{m) 


(1.3) 


The rest of this paper is organized as follows: In section 2 we introduce two different models 
called “inhnite-range” and “nearest-neighbour” models which are different in the range of their 
interactions. In section 3 we report the results of calculating the reduced density matrix of generic 
number of helds and compute entanglement and Renyi entropies of these two models. In section 4 
we investigate different features of these models probing them by entanglement measures including 
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entanglement inequalities and n-partite information. In the discussion section we will give some 
comments about the holographic dual of such a construction and also a field theoretic counterpart 
for black-hole evaporation process. In appendix A we explain some details related to the calculation 
of the reduced density matrix of our models. 

2 Kinetic Mixing Gaussian Models 

In this paper we are interested in Gaussian models as the simplest examples of interacting field 
theories which are analytically tractable. The most general wave functional for such models is given 
by [19] 

^[4>i] = AT exp I - ^ J ^ (t)i{x)Gij{x,y)^j{y'^, (2.1) 

where AA is a normalization constant and Gjj(rE,y)’s are complex valued functions which are sym¬ 
metric on i,j indices and also on the variables x and y. The corresponding (total) density matrix 
is constructed as ptot.i'/'i, = 'k*[())(]'I'[(/)j]. One can define a generic reduced density matrix by 
integrating out (without loss of generality) the first m number of the fields on the whole space-time 
as 

P(m)Wm+lAm+l,--- A'nAn] = j V(j)i ■ ■ ■ (2-2) 

where (c/)'^, • • • , cj)'^) is identified with {(pi, • • • , (pm) in the integrand. 

Since we are interested in analytically tractable simple models, in what follows we have chosen 
the same value of coupling constant between our mutually interacting field theories which means 
all off-diagonal non-vanishing elements of Gij take the same value. We are mainly interested in 
two models that we define in the following subsections. In the first model, any (pi interacts with all 
other fields (pj with {i ^ j). This model is called infinite-range model.^ In our second model any 
field (pi interacts only with its nearest neighbours which are (pi±i. We consider this model with a 
periodic boundary condition in the field space and call it the nearest-neighbour model. See Fig.l 
for a geometric realization of these models in the field space. 

Since we are interested in Gaussian models, in both of our models we consider kinetic mixing 
terms as the interaction between the free scalar helds, thus we are always dealing with marginal 
couplings. Note that both of these models in the special case where the total number of fields is 
two {N = 2) reduce to the massless interaction model in [13]. 

®This terminology is borrowed from the literature of statistical physics where e.g. an Ising model which all sites 
interact with each other are called infinite-range interacting Ising model. Here “range” refers to the field space 
rather than the real space. We thank Ali Naji for introducing us with this terminology. We also thank Julien 
Vidal for bringing our attention to related references [20], where aspects of a closely related model, knows as the 
Lipkin-Meshkov-Click (LMG) model, has been studied previously. 
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Figure 1: Schematic plots of infinite-range model (right) and nearest-neighbour model (left) in the 
held space. The blue lines connecting different helds represents the existence of an interaction term 
between the corresponding helds. The inhnite-range model clearly has much more interactions than 
the nearest-neighbour model. 


2.1 Infinite-Range Model 

The inhnite-range model is dehned by the following action 



where all ijij’s are interacting mutually with the same coupling constant A. The wave functional of 
this model is given by Eq.(2.1) where 


G{x,y) 


W{x,y) 
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(2.4) 


and W{x,y) = V ^ . We have briehy explained some details of this model in appendix 

A. 

One can easily show that this model can be diagonalized with the following eigenvalues 

= , a = 1 , 2 ,--- ,iV-l , = l + (iV-l)^, ( 2 . 5 ) 

and after the corresponding orthogonal transformation one can rewrite this model in terms of new 
(primed) degrees of freedom 

1 r ^ 

^=2 ( 2 - 6 ) 

i=i 
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It is an easily task to check that the positivity of the Hamiltonian restricts the value of A to the 
following window 


2 

N - 1 


< A < 2, 


(2.7) 


which we will consider in what follows as the range where this model is well-defined. This model is 
shown schematically in the field space in the right part of Fig.l. 


2.2 Nearest-Neighbour Model 

The nearest-neighbour model is defined with the following action 


S 


N 

/ d^x ^ + A ^ d^(l)id^(t)j 

i=l (i.i) 


( 2 . 8 ) 


where {i,j) means that the summation runs over two neighbours of each which are (j)i±i. Because 
of symmetry considerations we impose a periodic boundary condition (in the field space) such that 
the nearest neighbours of (pi are (p 2 and (pM- It is obvious that the number of interactions in this 
model is much less than the infinite-range model. The wave functional of this model is also given 
by Eq.(2.1) where 


G{x,y) 


W{x,y) 
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(2.9) 


and again W{x,y) = V~^ (see appendix A). 

One can easily show that the nearest-neighbour model can also be diagonalized and expressed 
in terms of new (primed) free fields just as Eq.(2.6). The eigenvalues of G for the case of = 2 is 


^1,2 = 1 T 


( 2 . 10 ) 


and for the case of N{> 2) is 


= 1 T A, ^2,3 = 1 =F A cos , AAr_ 2 ,Ar_l = 1 =F A COS 

Ai = l + X, A 2,3 = 1 - a cos , An- 1 , N = 1 - A cos 


N: even 
N: odd 


( 2 . 11 ) 


After performing the orthogonal transformation which leads to Eq.(2.6), one can compute the 
Hamiltonian of this model and show that the positivity of the Hamiltonian restricts the value of A 
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to the following windows 


-2 < A < 2 N-.2 

< — 1 < A < 1 N\ even (2.12) 

—1 < A < (cos ^ N-. odd 

In what follows we consider the above range for the coupling constant A where this model is well- 
defined. The schematic plot of the nearest-neighbour model is given in the left part of Fig.l. 


3 Entanglement and Renyi Entropies 

In this section we report the results of computing the reduced density matrix and hence entanglement 
and Renyi entropies in our models using replica trick. Here we skip the details of the messy 
calculations leading to Tr[/9”^^], and we just present the final results. The interested reader may 
find some details about the essential steps of the computations in appendix A. 


3.1 Infinite-Range Model 

Considering the infinite-range model one can show that using the definition of the reduced density 
matrix given in Eq.(2.2), together with the standard method of replication one can calculate 
Tr[/9”^)] which leads to (see appendix A for details) 


Tr 



=-^nn 


1 -|- f{m, N) cos 





i 


(1 - 

(i+a" 


(3.1) 


where 


f(m N) = _ m-m)Y{m) _ 


Y{m) 


1 /Ay 2m 

4V2y ' 2 +{m-l)X' 


(3.2) 


and is 

/(m,Ai) = Y^^. (3.3) 

Note that the normalization constant M plays no role in entanglement and Renyi entropies thus we 
will ignore it in what follows. Also note that in what follows we drop the index i of which regards 
to the discretized real space since all ^/s have the same value denoted by 

Since the m traced out fields together with the rest {N — m) fields build up the whole system 
(the total density matrix corresponds to a pure state), one would expect the above expression to be 
invariant under m —>■ (A^ — m) which is manifest in the expression of f{m,N). 

Now we are equipped with everything needed to apply the definitions given in Eq.(1.3) for 
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entanglement and Renyi entropies. One can read the entropies as 


s'") ^ x; »'“>({) = »(”)({) ^ 1 , S ^ ^ ,5(0 = »(£) ^ 1 , 

* * k^O 


) _ nln(l -g) -ln(l 


= 


1 — n 


> 40 = 


-ln(l -0 - 




1-e 


ln(^ 


(3.4) 


where the infinite sum is UV divergent. In order to regularize these expressions we use a smooth 
momentum cut-off, i.e., If we consider the {d — l)-dimensional spatial manifold to be a 

{d — l)-torus with size L, the infinite sum simplifies to 


d-i 


k^o \Mo 


Cd,d-1 



+ Cd,d-2 



+ • • • + Cdfi, 


(3.5) 


where Cj’s are constants that only depends on d. All the terms of the resultant entanglement 
entropy are divergent and depend on the UV cut-off e except the last one which is a universal 
term. To investigate the physical features of this model^ in the following sections we will consider 
this universal term which is proportional to Cdfi- Also note that according to Eq.(3.4) the whole 
A-dependence of entropies in this model is carried by (A) and s(A). See Fig. 2 where the universal 
part of entanglement entropy of this model is plotted for different values of m and N. 

Since in this paper we are dealing with entanglement in the field space, in what follows by 
entanglement and Renyi entropies we mean the “density” of these quantities which is defined as the 
entanglement and Renyi entropies in units of the infinite volume factor Eq.(3.5). This is obviously 
also true for the case of other entanglement measures which we define in the following including 
mutual and tripartite information. Thus here we have constructed the entanglement measures to be 
finite by definition. This is different from what happens in the case of spatial entanglement entropy. 
In that case some entanglement measures e.g. mutual information is defined by the whole expression 
of entanglement entropy which includes an area divergence but the divergent terms cancel out as 
long as the entangling regions do not have an intersection.® 


3.2 Nearest-Neighbour Model 

Next we consider the nearest-neighbour model which again by using the definition of the reduced 
density matrix given in Eq.(2.2) together with the standard method of replication we calculate 

“^This argument is also valid for the nearest-neighbour model. 

^Although in the case of spatial entanglement entropy it is well known that tripartite information and in general 
n-partite information with n > 2 are UV finite quantities even when the entangling regions share boundaries, this 
seems not to be generally correct as there is a counter example with corner shape entangling regions which have a 
single point as a common boundary [21]. 







s- s, 




Figure 2: Entanglement entropy of the infinite-range model as a function of coupling A. 
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Figure 3: Entanglement entropy for nearest-neighbour model as a function of coupling A. 


Tr for m neighbour fields out of N ones which leads to the following results 




i r=l 


1 + 


/ 27rr\ 

_COS I _ 1 

2Y-{m)g-{N — m — 1) — g-{N — m + 1) \ n J 


2Y{m)g-{N — m — 1) 


n 

S=1 


1 + 


2Yj^{m)gj^{N — m — 1) 


f2'Ks\ 


cos - 


2Y^{m)g+{N — m — 1) + g+{N — m + 1) \n J 


(3.6) 


where 


K 

2 


9 ±{ n )=n 


S=1 


), , (d±{N) + 2s 

1 — A cos ---TT 


V N 


Y{m) = -4 ( - 


(-0 


m+l 


Z{m + iy 


Y±{m) = Y(m) zb Irf(m), 


2 / 2iV -b 1 =F 1 


TT - 1, 


d±{N) = sin 


dim) 4 2 J ' Z(m + 1)’ 


m—1 


z(m)= n [ 


r=l 


1 — A COS ( —TT 
m 


(3.7) 


We define ^+(0) = ^ for consistency with the infinite-range model in the case of = 2.® Again 
®Note that the inhnite-range and nearest-neighbour models are the same for the case of = 2 and N = 3. 
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Figure 4: Renyi entropies for the infinite-range model (left) and nearest neighbour model (right). 


M is irrelevant to the calculation of entanglement and Renyi entropies. The result Eq.(3.6) is valid 
for 1 < m < N — 1. For the case of m = iV — 1 one should use the following expression 


n 

Tr [/3’"(m = iV - 1, iV)] = AT 

r=l 


1 + 


2Yd{m) 

— -cos 

2Yd{m) + 1 



Ydim) = 


X^ g+iN-2) 
16 g+{N) 


(3.8) 

which of course is equal to the result of m = 1 from Eq.(3.6) as expected. R is not hard to show 
that one can sum up the results of Eq.(3.6) and Eq.(3.8) in a single formula as 


Tr [p’"(m, N)] = (Tr [p’^(m, N)] Tr [p^{N - m, A^)])^ 


(3.9) 


which is valid for 1 < m < A^. The advantage of using this more compact formula is two-fold: it is 
no longer a piecewise formula and also the m ^ N — m symmetry becomes manifest in this form. 
Mathematically there is no difference between using Eq.(3.6) together with Eq.(3.8), or Eq.(3.9). 
In what follows we will continue with the first choice. 

The expressions for the entanglement and Renyi entropies are similar to the infinite-range model 
given in Eq.(3.4), and we just have to replace s^"')(^) and s(^) with s*^"')(^+) -|- s("'^(^_) and s(C+) + 
s(^_) respectively where are solutions of 

^_ 2Y±{m)g±{N - m - 1) _ 

1 -b 2Y±{m)g±{N — m — 1) ± ^^^(A^ — m -|- 1) 

Eor the case of m = A^ — 1 we consider s^”'^(^) and s(f) where ^ is defined as ^ 
note that as we have mentioned before, the structure of the regularization is independent of the 
interaction terms, thus in this model it exactly obeys the same structure of the previous model 
given in Eq.(3.5). 
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4 Aspects of Field Space Entanglement 


In this section we investigate some important features of these models based on the entanglement 
measures computed in the previous section. First we discuss about some features of entanglement 
and Renyi entropies of these two models. Next we study some physical constraints on entanglement 
measures which are known as entanglement inequalities. We also study n-partite information for 
certain values of n, and entanglement negativity as two other entanglement probes in our models. 
This analysis may be helpful to gain a more physical intuition about the structure of entanglement in 
these models and perhaps more generally some generic physical features of field space entanglement. 

4.1 Infinite-Range Versus Nearest-Neighbour Model 

In this subsection we are going to compare the infinite-range and the nearest-neighbour models using 
some graphical analysis. Previously in Fig. 2 and Fig. 3 we have plotted the entanglement entropy 
of these two models as a function of the coupling constant A. Note that the Hamiltonian positivity 
condition for these models which was given in Eq.(2.7) and Eq.(2.12), results in a iV-dependence for 
the valid range of coupling A. This has caused some asymmetries in the entanglement and Renyi 
entropies under A —?> — A. Also note that in the case of A = 0, since the vacuum state of the these 
models reduces to a direct product state, there is no entanglement between the specihed degrees of 
freedom in these models. Eig.4 shows the Renyi entropy for these models as a function of coupling 
A for various Renyi indices n. These plots clearly show that Sj jj is a decreasing function of n as 
expected. 

In Eig.5 we have demonstrated the m-dependence of the EE in these two models for three 
different values of A. Considering the coupling constant A, the domain of validity of the infinite- 
range model is wider than the nearest-neighbour model (compare Eq.(2.7) and Eq.(2.12)). As the 
value of A starts increasing from A = 0, for N > 3 which the distinction between these two models 
makes sense, the nearest-neighbour model reaches its maximum value of coupling constant, which 
we call Aji^ax) before the infinite-range one < -^max)- Since as A —Amax the maximum 

value of the corresponding EE diverges, the value of the EE for the nearest-neighbour model starts 
to grow much faster than the infinite-range one as A ^ -^max- Therefore there always exists a 
A*(< Aj,^ax < -^max) where the value of the EE of the nearest-neighbour model touches the value of 
that of the infinite-range one and gets larger values for A > A*. 

It is also interesting to study Renyi entropy as a function of Renyi index n. This is done in Fig. 6 
where we have plotted the Renyi entropy (normalized by entanglement entropy) in our models for 
various parameter values as a function of re. In this figure the dashed black curve corresponds to the 
value of entanglement entropy which coincides at re = 1 with Renyi entropy at arbitrary coupling 
A. There exists two other interesting limits of Renyi entropy corresponding to re —?• 0 and re —)■ oo. 
In the re —>■ 0 limit, one can easily check that Renyi entropy by definition, Eq.(1.3), reduces to the 
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Figure 5: Entanglement entropy of infinite-range (blue) and nearest-neighbour (orange) models for 
different values of coupling. From left to right: A = 0.8, A = 0.99 = A*, A = 0.999 and we have set 
N = 100. 

Hartley entropy 


= lim = logP, (4.1) 

n->-0 

where P is the dimension of the image of the reduced density matrix. Since in our models V is 
infinite, as it can be seen in Fig. 6, the Hartley entropy is divergent in this case. On the other hand 
in n ^ oo limit one finds the min-entropy 

hm ,SW = -logA^ax, (4.2) 

n—>-cxD 

where Amax is the largest eigenvalue of the reduced density matrix. In this case according to Fig. 6 
the Renyi entropy saturates to a constant value which depends on the value of the coupling A, as 
expected. Also note that in all cases the Renyi entropy is a decreasing function of the Renyi index 

n. 


4.2 Entanglement Ineqnalities 

In a general quantum-mechanical system or quantum field theory, entanglement entropy (and other 
measures of quantum entanglement) are proved to satisfy various inequalities. As a first example 
of such inequalities, we consider those dealing with Renyi entropy which was defined in Eq.(1.3). 
Renyi entropies must satisfy a variety of different inequalities such as [22] 

As(»)<0. |.(,„_1)S(«))>0, 

> 0. |A((„-1)S<»))<0. (4.3) 

As we mentioned before, the first inequality which shows Renyi entropy is a decreasing function of 
Renyi index n is satisfied in our models (see Fig. 6). It is a straight forward exercise to show that 
the other three inequalities are also satisfied in both of our models. 

In what follows in this subsection we consider other important inequalities which is expected to 
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Figure 6: Renyi entropy of infinite-range (left) and nearest-neighbour (right) models as a function 
of Renyi index n for N = 8 and m = 4. The dashed black curve corresponds to the value of 
entanglement entropy. 

be satisfied generally, based on the classification given in [23]: 

1) Sa> 0 (positivity of EE) 

This is a trivial property which we have checked it for different points in the parameter space 
of our models in the previous section (see Fig. 2 and Fig. 3). 

2) Sa + Sb > Saub (Subadditivity) 

This property can be rephrased in terms of the positivity of mutual information (MI) which is 
defined as^ 

I{A,B) = Sa + Sb-Saub. (4.4) 

MI is a quantity which measures the amount of shared information between A and B. While 
dealing with SEE, where A and B correspond to spatial subregions, MI is a UV finite measure of 
entanglement in contrast to EE. Clearly the subadditivity property implies the positivity of MI, i.e., 
I{A,B) > 0. Using the definition of Renyi entropy, one can also define mutual Renyi information 
(MRI) from the corresponding Renyi entropies as 

B) = 4 ") + SP - 4 "^^. (4.5) 

While dealing with SEE, it is known that MRI does not have a definite sign. It might be interesting 
to verify this property in the case of FSEE. 

Since we are dealing with FSEE, the Hilbert space decomposition we chose implied I{mi, m 2 ) = 
Smi + Sm 2 — <S'mi+m 2 ) where Srm is the FSEE for the case which we have integrated out {N — mi) 
fields (similarly for MRI). We have plotted MI and MRI for both the infinite-range and the nearest- 
neighbour models in Fig. 7 and Fig. 8 where we have considered the A and m-dependence of these 
quantities. In both of these figures, the blue curve corresponds to the case of MI, and other curves 
correspond to higher Renyi indices, i.e. MRI. MI is shown to be always positive in our models. It is 
^Note that the definition of MI does not restrict subsystems A and B to be complements. 
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Figure 7: Mutual Renyi information of infinite-range model. Left: MRI as a function of coupling A 
for N = 4, mi = l,m 2 = 2 and different value for Renyi index. Right: MRI as a function of m for 
A = 0.9 and N = 50, mi = m, m 2 = 50 — m with the same value of Renyi index. 
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Figure 8: Mutual Renyi information of nearest-neighbour model. Left: MRI as a function of coupling 
A for N = 4, mi = l,m 2 = 2 and different value for Renyi index. Right: MRI as a function of m 
for A = 0.9 and N = 50, mi = m, m 2 = 50 — m with the same value of Renyi index. 
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Figure 9: Intrinsic entropy of infinite-range and nearest-neighbour models as a function of coupling 
A for = 6, mi = 1 and different values of m 2 . 


worth to note that we could not find any region in the parameter space of the infinite-range model 
where the MRI admits negative values. The typical behavior of this quantity is similar to what is 
shown in Fig. 7 for specific values of the parameters. In the nearest-neighbour model the MRI have 
both positive and negative values as shown in Fig. 8. Note that while we deal with mi and m 2 which 
are complements, we expect the MRI to be symmetric with respect to half of the whole number of 
fields denoted by N (see the right plots in Fig.7 and Fig.8). 

3) Sa < Saub + Sb (Araki-Lieb inequality) 

This property which is also called the triangle inequality implies the positivity of the intrinsic 
entropy which is defined as 

J{B,A)=Saub + Sb-Sa , JiB,A)>0. (4.6) 

Some specific examples of this inequality in our models are depicted in Fig. 9. 

4) Saubuc + Sb < Saub + SbuC , Sa + Sc < Saub + SbuC (Strong subadditivity) 

Both of these inequalities are called strong subadditivity (SSA) and must hold in any quantum 
system. These inequalities physically mean that mutual information and intrinsic entropy must 
increases under inclusion. These inequalities hold in our models as we have plotted explicit examples 
of them in both of our models in Fig. 10. 

5) Sa + Sb + Sc + SaubuC < Saub + SauC + Sbuc (Monogamy of mutual information (MMI)) 

In spite of previously mentioned inequalities, which are general properties of entanglement mea¬ 
sures in any quantum system, MMI does not necessarily hold in any quantum system and thus it 
is not considered as feature of entanglement entropy. Again this inequality can be rephrased as the 
negativity of tripartite information, i.e. I^^\A,B,C) < 0, which is defined as 

(A, B, C) = Sa + Sb + Sc — Saub — Sauc — Sbuc + Saubuc (4.7) 

= I{A, B) + I{A, C) - I{A, B U C). 
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Figure 10: SSA inequalities in infinite-range model (left) and nearest-neighbour model (right) as a 
function of coupling A for A" = 8 , mi = 1, m 2 = 2, m 3 = 3. 



Figure 11: Tripartite information for infinite-range and nearest-neighbour models as a function of 
coupling A for N = 10, mi = l,m 2 = 3 and different values of m 3 . Note that 1^1 is always positive 
and in the latter case saturates to zero. 

Generally in quantum mechanics or even in QFTs, depending on how the Hilbert space is parti¬ 
tioned, can be positive, negative or zero. In Fig. 11 we have plotted for both of our models 
corresponding to different partitioning of the field space. As is shown in Fig.ll, this inequality 
does not hold in both of our models and more interestingly the tripartite information is always 
non-negative in these models. It is also interesting to note that in the case of mi -|- m 2 -|- m 3 = N 
the tripartite information becomes zero. According to second equality of (4.7) this is a reminiscent 
of models which exhibit extensive mutual information property [24]. 

4.3 n-partite Information 

In the context of quantum information theory, partitioning the system into re-parts, a new quantity 
known as re-partite information® is defined as [25] 

n n n 

^ X] ~ X] ^ SAi\jAj\jAk — .+ (-1)"'5'aiUA2U-UA„- (4.8) 

i=l i<j i<j<k 

®This n has nothing to do with the index of Renyi entropy. 
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Figure 12: 4-partite information for model I and II as a function of coupling A for N = 10, mi = 
m 2 = m 3 = 1 and different values of m 4 . 

It is obvious that according to the above formula the definition of 1-partite and 2-partite information 
reduce to EE and MI respectively. Also note that the n-partite information for n > 1 is a UV finite 
quantity. Actually a finite measure for quantum entanglement between subsystems of a larger 
system is not unique (e.g. another choice known as multi-partite information is defined in [27]). 
The reason why we use the above definition for n-partite information Eq.(4.8) for such a quantity 
is due to its property which reduces to the definition of tripartite information Eq.(4.7) in the case 
of n = 3 (while e.g. multi-partite information does not have this property [27]). 

As we have mentioned before, MI is always non-negative, i.e., > 0, due to the subadditivity 

property of EE. Although the sign of tripartite information is not fixed generally, but as we have 
shown in the previous subsection it is always non-negative in both of our models. It is worth to 
note that in the case of CFTs which support a gravitational dual, the sign of tripartite information 
is hxed to be always negative. This general property restricts the holographic mutual information 
to be monogamous [25].® As an extension of this property, it is also shown in reference [28] that in 
a specific limit in the case of SEE, the holographic n-partite information has a definite sign: it is 
positive (negative) for even (odd) n. 

It would be interesting to investigate the sign of higher n-partite information in our models. In 
Fig. 12 and Fig. 13 we present the 4-patite and 5-partite information as a function of the coupling A 
which is surprisingly always positive. Also focusing on 5-partite information together with 3-partite 
information (see Fig. 11), one may conjecture that n-partite information is always vanishing for the 
case of odd n’s with complement partitioning of the system i.e. m, = N. 

4.4 Entanglement Negativity 

Entanglement negativity and its counterpart logarithmic negativity are useful measures of quantum 
entanglement even for mixed states [30]. It is known that the von-Neumann entropy for a mixed 
state, e.g. a thermal state, dominated by the classical correlations is not a useful measure for 

®It is also shown in reference [26] that the null energy condition is a necessary condition for the monogamy of 
holographic mutual information. 
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Figure 13: 5-partite information for model I and II as a function of coupling A for N = 10, mi = 
m 2 = m 3 = m 4 = 1 and different values of m^. Note that is always positive and in the latter 
case saturates to zero. 

quantum entanglement. MI also measures the total correlations (both quantum and classical) 
between two subsystems which just offers an upper bound [29]. It has been shown that negativity 
is an entanglement monotone (does not increase under any LOCC operations) and hence a proper 
measure for quantum entanglement [31]. To give a more concrete but nevertheless simple definition 
of this quantity one may consider a tripartite system in a pure state with a complement partitioning, 
i.e., M = Ai D A 2 U A 3 . In this case the reduced density matrix corresponding to union of two 
subsystems is described by a mixed state p = paiuA 2 - Entanglement negativity and logarithmic 
negativity are defined as 

AA(p)^^^y^, T(p)=log[|/^||, (4.9) 

where [|p^^|| denotes the trace norm of the partial transpose of p. With the above definition the 
logarithmic negativity measures how much the eigenvalues of \\p'^^\\ are negative. 

Although computing these quantities in general is not an easy task, the authors of [34] have 
introduced a replica approach to obtain the logarithmic negativity in the ground state of 2d CFTs. 
They also show that for a pure state and bipartite system where H = Hi ^ H 2 , this quantity is 
given by Renyi entropy with n = 1/2, i.e., 

T(P2) = 21ogTVp^/^ (4.10) 

We focus on this definition in order to study the logarithmic negativity in our models. We postpone 
further investigations based on computing Eq.(14) for future works. In Fig. 14 we have plotted 
logarithmic negativity as a function of coupling A for different partitions of the Hilbert space. 
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Figure 14: Logarithmic negativity for infinite-range and nearest-neighbour models as a function of 
coupling A for = 10 and different values of m. 

5 Conclusions and Discussions 

In this paper we have considered a less studied type of entanglement which is known as field space 
entanglement. This type of entanglement corresponds to a Hilbert space decomposition in the field 
space of a quantum field theory. As a simple laboratory to study field space entanglement, we 
have considered a theory with a generic N number of free scalar fields, we added kinetic mixing 
terms (in terms of two specific models) which generates entanglement between these scalar fields. 
We traced out a generic m number of these fields and worked out the entanglement and Renyi 
entropies between m and (N — m) number of these scalar fields. The result of these entropies is 
UV-divergent which scales with the (spatial) volume of the theory as expected. Similar to the case 
of spatial entanglement entropy, there is a universal term, i.e. a UV cut-off independent term which 
we argue to carry some information about the theory. Beside the entanglement and Renyi entropies, 
we also constructed other well known entanglement measures such as mutual information, intrinsic 
entropy and n-partite information to further investigate features of field space entanglement. We 
have shown that this type of entanglement in our models satisfy most of the known general features 
of entanglement measures including Renyi entropy inequalities, strong subadditivity and Araki-Lieb 
inequality. We have also studied the monogamy of mutual information which has a definite sign 
(positive) for tripartite, 4-partite, and 5-partite information in our models. 

There are several directions which one can follow to further investigate our models and the 
notion of field space entanglement using this laboratory. We leave further investigations of these 
models, including the recently proposed entanglement inequalities (see [32]), to future works and 
in the following of this section we discuss a few words about the holographic picture of field space 
entanglement entropy and also offer a different viewpoint to this family of field theories which we 
have considered. 

Holographic Picture of FSEE 

In order to gain some information about the possible gravity picture of such an analysis, as the first 
step we consider some well known features of field theories which support holographic duals: the 
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monogamy condition for holographic mutual information and its implication on the dual field theory. 
As we mentioned in the previous section the tripartite information in both of our models is always 
positive and the monogamy constraint does not hold. Actually this behavior is in contrast to the 
holographic result which shows that the holographic mutual information is always monogamous [25]. 
So in this sense it seems that our models do not have a well defined holographic description. It is 
important to mention that it is not clear that whether this constraint must hold for any type of 
EE or it is just a feature of SEE. In the following for a while we forget about this comment on the 
relation between monogamy of mutual information and the existence of a holographic dual. 

The authors of reference [13] have proposed a naive holographic picture for the entanglement 
entropy between two CETs which might be related to our models in the case of = 2. In this 
proposal the factorization of the Hilbert space in the field space was related to partitioning the 
compact part of the AdSs x S® geometry by introducing a dA surface which partitions the S® sphere 
into two parts and wraps the boundary of AdSs. The minimal surface anchoring the corresponding 
boundary on a certain UV cut-off surface was proposed to give the entanglement entropy between 
two interacting subsectors of the whole CFT 4 (which is dual to the AdSsxS^ geometry). Although 
there are some substantive comments about the relation between this holographic picture and FSEE 
(see [ 6 ] and also [14]), the holographic dual of our models in this picture is straightforward. One 
may partition the S® sphere to N parts and the corresponding entanglement entropy is proportional 
to the volume of different portions. For example if we consider the mutual information between 
two set of fields, the sphere is divided into three parts and different terms contributing in the 
expression of mutual information are proportional to the volume of the corresponding part of the 
sphere. 

There is another geometrical picture introduced in reference [14] which offers a geometrical 
interpretation for the entanglement between two SU(m) and SU(A^ — m) CETs again as subsectors 
of the dual CFT 4 . This picture is based on the interpretation of minimal surfaces in the more 
general supergravity Coulomb branch geometry rather than AdSsxS® as entanglement entropies. 
Here the level sets of the scale factor multiplying the Minkowski part of the solution is interpreted 
as the UV cut-off of the CETs living on separated stacks of D3-branes. There are two family of level 
sets: disconnected level sets which are consisted of two separated surfaces surrounding each brane 
stack, and connected ones which are single surfaces surrounding both brane stacks. Correspondingly 
there are two family of minimal surfaces, those which start and end on the connected level sets and 
those which start and end on the disconnected level sets. Those surfaces which start and end on 
the connected level sets are interpreted as a measure for the entanglement between two CETs living 
on the brane stacks which is generated by means of the stretched modes between these stacks. The 
minimal surfaces starting and ending on a part of the disconnected level set around, say stack 1 , 
are interpreted as a measure for the entanglement between a part of CFTi and CFT 2 living on 
the other stack together with the entanglement between two parts of CFTi. For more details see 
reference [14]. 

One can naively generalize this picture to be appropriate for interpreting mutual information 
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between any two of three SU(mi) and SU(m 2 ) and SU(A'' — mi — m 2 ) CFTs by considering three 
stacks of D3-branes. In this case the number of connected and disconnected level sets increase. 
There are four types of disconnected level sets: a single one composed of three parts and those 
which are composed of two parts, one surrounding two stacks and the other surrounding a single 
stack. Although this configuration for three stacks is too complicated to calculate, there are several 
minimal surfaces which could be interpreted as a direct generalization of what was discussed in the 
previous paragraph. One can in principle even generalize this picture for arbitrary N and interpret 
the corresponding minimal surfaces as in the case of = 2 as a possible holographic picture of our 
models. 

On the other hand it is recently argued in reference [6] that it is not possible to give a precise 
geometrical realization for FSEE in a holographic dual and all which is discussed in the above two 
scenarios is rather related to entanglement in the space of the global symmetry of the CFTs which 
is in no way essential to define FSEE. Although the author has offered some arguments to give an 
effective realization to such a case in terms of IR CFTs as dual field theories for internal throats in 
the Coulomb branch supergravity solution of separated D3-branes, the geometrical interpretation 
for FSEE seems to still be an open problem. 

Now lets forget about different scenarios as candidates for the holographic picture of FSEE. One 
may focus on the A^-dependence of the entanglement entropy in the infinite-range model to give a 
concrete expectation for a possible reliable holographic dual.^*^ To avoid unnecessary complications, 
we consider the entanglement entropy in the leading order of A 


S{m) = — m{N — m) 
32 


1 - log 


X‘^m{N — m) 


32 


+ o(a3), 


which for the special case of m = ^ gives 


S{m) = 


A^Ai^ r 


128 


1 - log 


128 


+ 0(A3), 


(5.1) 


(5.2) 


which is expected to be explained by any holographic dual. One can work out the corresponding 
expressions for the nearest-neighbour model. 

Beside this check, the large N behavior of these models seems to have interesting features in 
the field space. In this limit the infinite-range model seems to behave as a non-local theory in the 
field space while the nearest-neighbor model resembles a local theory.It would be interesting to 
investigate this property more precisely and study its implications specifically on entanglement and 
Renyi entropies. 

thank Shahin Sheikh-Jabbari because of his valuable comment about the A^’-dependence of field space entan¬ 
glement entropy which was insightful for us to clarify the structure of our analysis. 

'^^We thank Shahin Sheikh-Jabbri for drawing our attention to this interesting point. 
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Figure 15: The “information” for the infinte-range model for N = 100 and A = 0.9. 

A Model for Black-hole Radiation 

A field theory which consists of a number of interacting fields could be a field theoretic counterpart 
of Page’s model for black-hole evaporation process [33].^^ A first and simple clue for this argument 
is the symmetric behaviour of the entanglement entropy around m = ^ (see Fig. 5 were we have 
plotted this behavior for both of our models) and one may compare it with the entanglement (or 
information) evolution during the black-hole evaporation. 

In reference [33] the author has considered two subsystems with Hilbert space dimensions m 
and n respectively such that the total Hilbert space with dimension m x n is in a pure state. 
He has shown that the entanglement entropy between these two subsystems is symmetric as a 
function of the thermodynamical entropy which is defined by logm. Another important result of 
such a consideration is that the deviation of the entanglement entropy from its maximum value 
(the thermodynamical entropy), which is defined as “information”, remains almost zero until the 
entanglement entropy reaches its maximum value. 

We demonstrate the entanglement entropy (see Fig. 5) and “information” (see Fig. 15) as a func¬ 
tion of m. The information is defines as I = m — S. Our argument for considering such a definition 
for information in this case is as follows; In our model where the total Hilbert space includes N 
fields, the subsystems (I) and (H) have m and (N — m) fields respectively and the thermodynamical 
entropy is an extensive quantity. To see this consider the Hilbert space for the first subsystem which 
is 'H(i) = 'Hi ®'H 2 ® so if we denote the dimension of the Hilbert space for a single field by 

T), then the dimension of "^( 7 ) becomes D™. So in our case the themodynamical entropy becomes 
log = m log D and we expect that in the definition of information one must replace log m with 
m. 

Note that in Fig. 15 which we have plotted the information I, it is non-zero even in the early 
stages of evolution (m ~ 1), in contrast with what was previously found in [33]. 

thank Mohsen Alishahiha for bringing our attention to this point. 
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A Calculation of Reduced Density Matrix 

In this section we explain some details of the calculation of our master formula, which is the trace of 
the reduced density matrix of both of our models reported in Eq.(3.1) and Eq.(3.6). Here we explain 
the logical steps with general formulas as the key points leading to these results. The remaining part, 
although is some how messy, it is of course straightforward if one follows the procedure discussed in 
this section. The starting point is the wave functional for Gaussian models introduced in Eq.(2.1). 
We explain the general formalism while explaining the infinite-range model in subsection A.l, and 
turn to the nearest-neighbour model in subsection A.2. 


A.l Infinite-Range Model 

As we have mentioned in Sec. 2, the total density matrix of these models is generally defined as 

• • • ,(I)'n^4>n]=^*[4>'i,4>2,--- - ■ ■ ,4>n], (A.l) 


where T[{</>}] is the Gaussian wave functional introduced in Eq.(2.1). In order to define the reduced 
density matrix for the simplest case, i.e. m = 1, we identify cpi and (/>( and integrate over it on the 
whole space 

P(N-102,<i>2,4>3i4>3r ■ ■ ,(Pn,4>n] = j ,<Pn]- (A.2) 

Implementing the explicit form of the Gaussian wave functional given in Eq.(2.1) and performing 
the integral, up to an irrelevant normalization constant the result is 


P(Ar-i)[(/>2,(?5'2,--- =exp| - ^ ^ J dx^ ^ 0x) 


Grj- 


GuGij 

01 




+ 0X) ( GF - -00 ) 0y) - Mx)GuGl0iy) - 0x)G*uG,j^,iy) 
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11 


(A.3) 
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where we have dropped the x and y dependence of Gjj’s in the above expression for simplicity and 
we do so in what follows. Note that in the above formula • = 2Re [•]. It is not a hard task to 
integrate out more than one field, say m number of fields which leads to the reduced density matrix 

P{N-m)Wni+lAm+ir-- A'nAn] = j V4>i ■ ■ ■ V(prn[4>1, ^^2 , ' ' ' A'n]^ [(pl, 4>2, ' ' ' ,4 >n]- (A.4) 

A similar procedure which leads to Eq.(A.3) can be performed to arrive at (via induction) 

P{N—m) exp \ ^ ^ / dx dy X 

I i,j=m^l 

'Ux)x\fcl^j{y) + cl>[{x)x\f*ct^'^iy) + Ux)Y^U'j{y) - 




(A.5) 


where 


Xij - 



Zim.) 


yim-l) 

YYmm 


yM _ y(").-l) 
b b 
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yim-l) 


_ y(m-l) -^(m-1) 
^ i,m—l ‘ 


(A.6) 

One can work out the generic reduced density matrix using the above recursion relations with initial 
values = Gij and = 0. Considering the inhnite-range model, using Eq.(A.7) together with 
Eq.(2.4), one can find the most general form of the reduced density matrix in terms of m, N and 
A which is the coupling constant between the scalar fields. Eor future use we rewrite the reduced 
density matrix as 


P(N—m)[4^m+li 4^'m+lj ' ' ' exp 



(pm +1 ( 3 ^) 


1 

2 


/ 


dx^-Hy<^-^ 


X 





M{m, N) 


(pm+l {y) 

(Pniy) 

\ (t>N{y) / 


(A.7) 
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Figure 16; Replica method for N = A and n = 3 for m = 1 (left) and m = 2 (right). The spatial 
directions of the field theory are perpendicular to the plane and the vertical lines correspond to the 
time direction. The numbers under each vertical line corresponds to i-th. field (pi. 
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(A.8) 


After the construction of the reduced density matrix, one can use the standard replica method 
[19,35-37] to construct the its n-th power in order to work out its trace. This step is basically the 
same for both of our models which is pictorially explained in Fig. 16 for m = 1 and m = 2 and 
= 4. The replica method here is exactly the same as the well-known procedure for 2d CFTs 
within the context of spatial entanglement (e.g. see [37]). The only difference is that here we cut 
along the whole spatial coordinates at r = 0 of those fields which we are not integrating out (see 
Fig. 16). 

What remains to do is to start from Eq.(A.7) and find the trace of the reduced density matrix 
for general Renyi index n for generic m and N. It is not a hard task, although messy, to see that 
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using replica method one can find 


Tr 
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(A.9) 


where the matrix A4 is a n{N — m) x n{N — m) square matrix and is defined 
blocks as 

M 

■M.m+2,m+l ■^m+2,m+2 M m+2,N I 

M= _ . . . > 

\ Mn ,m+l MN,m+2 ■ ■ ■ Mn,N / 

and the blocks M.rn,m' are n x n square matrices given by 
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(A.ll) 


If we calculate the determinate of M we are done. This would be a much simpler task if we 
consider the explicit values of Gjj’s for the inhnite-range model. To do so the key point is the 
existence of an orthogonal transformation which results in a diagonal model (free scalar helds) as 
was explained in section 2.1 and specifically in Eq.(2.6). In the diagonal basis the ground state 
wave functional up to a normalization constant becomes 


^[<(>' 1 , • • • A'n]= exp 
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2 
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W{x,y) 


Ai(l)'i{x)(j)'i{y) 


(A.12) 
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where Aj’s are given in Eq.(2.5) and W{x,y) is given by 


k 


(A.13) 


where V is the {d — l)-dimensional volume which the field theory is defined on. Since we have 
applied an orthogonal transformation between {(pi, • • • , ^at} and {cp'i, • • • , basis, the physical 
state is unaffected, i.e. 

^Wi,--- A'n\ = An], 

and we can rewrite the ground state in terms of {(pi, ■ ■ ■ , cpM^ basis as 


,(pN]=eyiv{-- j dx dy W{x,y) 
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Gij(pi{x)(pj{y) 

*j=i 


(A.14) 


where Gjj’s for this model are given by 
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(A.15) 


Using these explicit expressions and working out the trace of the reduced density matrix hrst for 
m = 1 and generic N, by induction one can easily find that 
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(A.16) 


Now we are done with the m = 1 case. Generalizing to m > 1 is not a hard task because of a 
simple structure in the reduced density matrix. Since the structure of the reduced density matrix 
only depends on [N — m) rather than m and N itself, we are almost done since we already have 
calculated m = 1 for generic N. Again by induction one can generalize the above result for general 
m which is 


Tr 
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where Y{m) (not to be confused with the Y elements of matrix A4) is defined as 
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(A.18) 


27 


















A.2 Nearest-Neighbour Model 


The logical steps for this model is the same as that we have discussed in the previous subsection. We 
may start from Eq.(A.4) for this model. In comparison with the infinite-range model, this model 
has much fewer symmetries which makes it harder to push this calculation as general as we did 
for the infinite-range model. Since we are interested in the case where the strength of interactions 
between interacting fields is equal, we will restrict our analysis for equal off diagonal values of Gij 
which we denote by Gij = G ioi i ^ j, and also Gu = G^- For such a case one can perform m 
number of Gaussian integrals to arrive at the general reduced density matrix in the form of Eq.(A.7) 
with 


M(m, N) 


/ V* ^ 

\Xm\^ 

G* 

0 

^ m 

\Ym? \ 


Gd-Xl 

0 

G 

|y |2 

1 m 1 

Y^ 

G* 

0 

Gd 

0 

0 

0 

0 

0 

0 

Gd 

0 

0 

^ m 


0 

0 

2 

\Xm\^ 

V 

Y^ 

0 

0 


Gd - X^ j 


where 


Xm=G 




Ym=G 


G 


m—l 


{-2y 


nT=i 


Zm — Z\ — 


G2 




(A.19) 


(A.20) 


and Z\ = Gd- Note that the above general form is correct for m > 1, for the case of m = 1 there is 
an extra factor of 2 in the denominator of all components represented in terms of Ym- The reader 
should note that these and Z^ functions are not to be confused with the functions with Y(m}, 
Yd(m) and Z(m) which appear in the final result as functions of the coupling which is given in 
Eq.(3.6) and Eq.(3.7). 

Now we can work out the counterpart of Eq.(A.9) in this model. Here the form of A4 is more 
complicated and is given as follows 


/M^ 


0 

0 


M^\ 



MP 

0 


0 

0 


Jvpd 



0 

0 

0 



Mpd 


\m^ 

0 


0 


M^J 
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where again the blocks M’’ are n x n square matrices given by 


fCd-Xl \x„ 


M^ = 


l^ml 


Gd-Xi 

\Xm\^ 


0 
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V \Xm? 

( n 0 

|y^|2 


M^ = 


0 

Gd-Xl 

\X^? 
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• • • lyr. 


\Xm\^ \ 

0 
0 


Gd - XI \X„ 


1^. 


Gd-XlJ 


lA 


0 \Yrr, 


Vi^;. 


y2 

iy I' 

I m I 
0 


y2 

m 


\Yrr, 


0 

0 

|y b 

m\ 

y2 

m 


(A.22) 


and = diag{G, • • ■ , G} and M'^<^ = diag{Grf, • • ■ , G^}. 

Now we are equipped with Tr for the nearest-neighbour model and what remains is 

to plug in the corresponding Gij which was given in Eq.(2.9) and work out the determinant of Ai 
given in Eq.(A.21). This step is of course more messy than the case of infinite-range model because 
of a technical subtlety. Here in contrast with the infinite-range model, when we increase m and 
N, the degree of the polynomials appearing in the expression of det[Af] also increases. The key 
point to bring these expressions back into control is to factor them in terms of their roots, which 
generally take the form of = cos [w{m, iV)7r] with different w{m, N) functions. Eollowing such 
a process will lead to Eq.(3.6). Note that the functions X and Y used here has nothing to do with 
the functions given in the hnal result Eq.(3.6). 
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